Cosmic Masers - from OH to Ho 
Proceedings lAU Symposium No. 287, 2012 
R. Booth, E. Humphries & W. Vlemmings, eds. 



@ 2012 International Astronomical Union 
DOI: 00.0000/XOOOOOOOOOOOOOOOX 



OH Masers and Supernova Remnants 

Mark Wardle and Korinne McDonnell 

Department of Physics & Astronomy and Research Centre for Astronomy, Astrophysics & 
Astrophotonics, Macquarie University, Sydney NSW 2109, AustraUa 
email: mark . wairdleQmq . edu . au , korinne . mcdonnellOmq . edu . au 

Abstract. OH(1720 MHz) masers are created by the interaction of supernova remnants with 
molecular clouds. These masers are pumped by collisions in warm, shocked molecular gas with 
OH column densities in the range lO'^*' -10^*^ cm^^. Excitation calculations suggest that inver- 
sion of the 6049 MHz OH lino may occur at the higher column densities that have been in- 
ferred from main-line absorption studies of supernova remnants with the Green Bank Telescope. 
OH(6049 MHz) masers have therefore been proposed as a complementary indicator of remnant- 
cloud interaction. 

This motivated searches for 6049 MHz mascr emission from supernova remnants using the 
Parkes 63 m and Effelsberg 100 m telescopes, and the Australia Telescope Compact Array. A total 
of forty-one remnants have been examined by one or more of these surveys, but without success. 
To check the accuracy of the OH column densities inferred from the single-dish observations 
we modelled OH absorption at 1667 MHz observed with the Very Large Array towards three 
supernova remnants, IC 443, W44 and 3C 391. The results arc mixed - the OH column is revised 
upwards in 1C443, downwards in 3C391, and is somewhat reduced in W44. We conclude that 
OH columns exceeding lO^'^ cm~^ are indeed present in some supernova remnants and so the 
lack of any detections is not explained by low OH column density. Wo discuss the possibility 
that non-local line overlap is responsible for suppressing the inversion of the 6049 MHz line. 

Keywords, masers, line: formation, molecular processes, shock waves, ISM: molecules, super- 
nova remnants, radio continuum: ISM, radio lines: ISM 



1. Introduction 

OH (1720 MHz) masers associated with supernova remnants are unambiguous indica- 
tors of interaction with adjacent molecular clouds (Frail, Goss & Slysh 1994). The relative 
ease of identifying supernova remnant - molecular cloud interaction sites relative to other 
methods, such as searching for molecular gas that is kinematically disturbed or chemi- 
cally processed (e.g. DeNoyer 1979a,b; Wootten 1981) motivated surveys of the galaxy 
for more examples (Frail et al 1996; Yusef-Zadeh et al 1996; Green et al 1997; Koralesky 
et al 1998; Yusef-Zadeh et al 1999a,b; Sjouwerman & Pihlstrom 2008), which have found 
that approximately 10% of galactic supernova remnants have 1720 MHz masers and are 
therefore interacting with molecular clouds. 

The masers are collisionally pumped in warm, dense molecular gas with densities and 
temperatures in the range 50-125 K, ~ 10^-10^cm~^ and OH column densities 10^^- 
10^^ cm~^ (Elitzur 1976; Lockett, Gauthier & Elitzur 1999). These conditions are not met 
in gas overrun by standard "J-type" shock waves which dissociate molecules: although 
molecules reform as the shocked gas cools (e.g. HoUenbach & McKee 1989), cooling 
is so rapid that the OH column that is produced in the relevant temperature range is 
too small (Lockett, Gauthier & Elitzur 1999). Instead the molecular shock must be "C- 
type" (MuUan 1971; Draine 1980), in which heating and compression of the gas within the 
shock front is achieved by magnetically-driven ion-neutral collisions. The low fractional 
ionisation of the molecular gas means that this process is relatively slow, yielding a broad 
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Figure 1. OH absorption from ATCA observations of SNR W28 (Green, Wardle & Lazendic 
2000), overlaid on the locations of bright OH(1720 MHz) maser spots (Claussen et al. 1997), 
and contours of CO J=l-0 and J=3-2 emission (Arikawa et al. 1999) tracing unshocked and 
shocked gas, respectively. (Adapted from Arikawa et al. 1999) 

shock front with an extended warm tail. However the requisite OH column is still higher 
than predicted by standard C-type shock models which rapidly convert OH into water 
(Draine, Roberge & Dalgarno 1983; Kaufman & Neufeld 1996). Instead an additional 
agent, such as X-rays or cosmic rays associated with the supernova remnant, is needed 
to dissociate water into OH as the gas cools behind the shock, (Wardle 1999). 

The clear association of 1720 MHz mascrs with shocked gas (see e.g. Fig. [ij, under- 
pinned by a theoretical production mechanism hinging on molecular shock waves, pro- 
vides a firm basis for assuming that the presence of these masers signals the existence of a 
shock wave driven by a supernova remnant into an adjacent cloud (Wardle & Yusef-Zadeh 
2002). This has been confirmed by follow-up studies of individual remnants. A good illus- 
tration of this is provided by the supernova remnant G349.7-I-0.2, about which little was 
known apart from its having a radio continuum shell and four OH(1720 MHz) masers. A 
search for CO (J=l-0) emission from an adjacent molecular cloud, H2 (v=l-0) emission 
at 2.12 /im from shocked molecular gas, and OH absorption associated with the masers 
found all three (Lazendic et al. 2010). The population of 1720 MHz-bearing supernova 
remnants show other hallmarks of interaction. For example "mixed-morphology" rem- 
nants are characterised a radio shell, or partial shell, filled with centrally-peaked emission 
in soft thermal x-rays (Long et al. 1991; Rho & Petre 1998) believed to be provided by 
the evaporation of dense interstellar gas in the interior of the remnants (White & Long 
1991; Cox et al. 1999). There is a strong statistical association between remnants hosting 
1720 MHz masers and mixed-morphology remnants (Yusef-Zadeh et al. 2003), and also 
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with those with associated gamma-ray sources (Hewitt et al. 2009), which are hkely pro- 
duced by interactions of cosmic-rays accelerated in the remnant with the neighbouring 
molecular gas. 

Apart from signalling locations where remnants are driving shock waves into molecular 
clouds, 1720 MHz masers are important probes of the interaction. The masers require 
significant columns of shocked molecular gas and good velocity coherence, and so tend 
to occur on the limbs of supernova remnant shells at velocities close to the systemic LSR 
velocity of the remnant, allowing distance estimates based on models of galactocentric 
rotation. Follow-up Zeeman observations at 1720 MHz have determined that the magnetic 
field in the post-shock gas is ~ 1-3 mG (Claussen et al 1997, 1999; Brogan et al 2000; 
Hoffman et al 2005). The post-shock magnetic field acts as the piston driving a C-type 
shock wave, and a field of this strength provides sufficient post-shock pressure needed 
to drive a shock at 25kms~^ into gas of density ~ lO'^cm"^, so these measurements 
provide an important sanity check on the basic picture. In addition the magnetic field 
measurements provides a good estimate of the pressure in the interior of the supernova 
remnant that is driving the shock into the molecular cloud. 

2. Searches for 6049 MHz masers 

While 1720 MHz masers are invaluable in pinpointing where supernova remnants in- 
teract with molecular clouds, their absence docs not imply that an interaction is not 
occurring. Instead, one or more of the physical conditions required to collisionally invert 
the 1720 MHz transition may simply be absent - a sufficient column of OH, or line-of- 
sight velocity coherence within the OH column. Additional indicators of such interactions 
would therefore be useful. 

The 1720 MHz inversion dies away for OH columns > 10^'' cm~^ due to photon trap- 
ping, but under similar physical conditions the analogous 6049 MHz transition in the first 
rotationally-excited state of OH becomes inverted (Pavlakis & Kylafis 2000; McDonnell, 
Wardle & Vaughan 2008). This prompted searches for masers at 6049 MHz associated 
with supernova remnants using Parkcs (McDonnell, Wardle & Vaughan 2008), Effcls- 
berg (Fish, Sjouwerman & Philstrom 2007) and the ATCA (McDonnell 2011) (see Table 
1). These searches targeted SNR with 1720 MHz masers, mixed-morphology, or other 
evidence of interaction (note that there is substantial overlap between the SNR in the 
Parkes & Effelsberg and Parkes & ATCA searches), but all were unsuccessful. 

Table 1. Searches for 6049 MHz masers 
Telescope A'snr beam a (mjy/beam) Reference 

Parkes 35 200" ~ 150 McDonnell, Wardle & Vaughan (2008) 

Effelsberg 14 130" ~ 10 Fish, Sjouwerman & Philstrom (2007) 

ATCA 17 20" X 9" ~ 22 McDonnell (2011) 



The non-detections can be converted into upper limits on the 6049 MHz maser optical 
depth (i.e. maser amplification factor) using estimates of the brightness temperature 
of the background continuum and the likely beam filling factor of any putative maser 
emission. In turn, an upper limit on the maser optical depth can be translated to an 
upper limit on OH column density through the OH excitation calculations. It turns 
out that the relatively poor sensitivity of the Parkes observations does not place strong 
constraints on the OH column; by contrast, the derived limits on OH column density 
implied by the lack of detected 6049 MHz masers in the Effelsberg and ATCA searches 
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are < 10^^'^ cm^-^ (McDonnell 2011). These upper limits are at odds with the OH columns 
exceeding ^ 10^^ cm~^ that were inferred in nine supernova remnants by modelling the of 
OH absorption in the 1612, 1665, 1667MHz lines observed with the Green Bank Telescope 
(Hewitt et al. 2006; Hewitt, Yusef-Zadeh & Wardle 2008). These nine remnants are listed 
in Table 2 - all but G349.7+0.2 were searched for 6049 MHz emission with Effelsberg 
and/or ATCA, and so one might have expected several detections. 

One potential source of error is the 7.4 arc minute beam of the GET observations. 
While a beam filling factor was included in the modelling by Hewitt et al. (2008) this does 
not account for different sub-structure in the background continuum and overlying OH 
absorption within the beam, and given that the inferred beam filling factors are typically 
< 10~^ this omission may be significant. To check whether the inferred OH columns are 
systematically overestimated we conducted 1667 MHz OH absorption observations of 3 
remnants (3C391, W44, and IC443) using the VLA in D array, yielding a beam size of 
^ 55" X 90". We followed the approach of Hewitt et al. (2008) in modelling the absorption 
at 1667 MHz in conjunction with 1720 MHz VLA data (the latter kindly provided by Jack 
Hewitt). 

Notably, the derived solid angles subtended by the OH absorption components ~ which 
are typically marginally resolved in the VLA observations - are in reasonable agreement 
with the GBT results. However, the comparison of the derived OH columns gives mixed 
results: the peak column density derived from the VLA observations for 3C391, W44 and 
IC443 are much less, a little less, and somewhat greater than the GBT-derived columns, 
respectively. 

Fig. [2] shows the 1667 MHz absorption against the continuum in IC 443. As expected, 
OH absorption surrounds the three known 1720 MHz OH masers, which are associated 
with molecular clumps clumps D and G and a new feature, labeled OHl. Focussing on 
clump G, which has the greatest column density, we derived model parameters quite 
similar to the GBT. This reflects that the OH emission is dominated by a single compact 
component and the background continuum is quite smooth. For 3C491 by contrast, the 
continuum shows strong structure on the scale of the dominant absorption in the NW 
of the remnant (see Fig. 3). In this case the GBT observations, which have a beam area 
of about the size of the image in Fig. |3j underestimate the strength of the background 
continuum underlying the absorption by averaging over the beam, resulting in an over- 
estimate of the optical depth in the 1667 MHz line. The derived OH column density is 
therefore too large. 

In summary, the VLA observations show that the GBT modelling of the does pretty 
well, but the poor resolution of single dish observations leads to significant overestimates 
of the OH column density when the background continuum and OH absorption have 
poorly correlated substructure within the beam. Nevertheless IC 443 and potentially 
some of the other remnants in Table 2 have OH column densities exceeding 10^^ cm~'^ 
and yet were not detected in the ATGA and Effelsberg 6049 MHz searches. 

3. Line overlap 

An alternative possibility is that the treatment of overlapping spectral lines in the OH 
excitation calculations is at fault. Line overlap in OH occurs even for velocity widths 
of order lkms~^ because of the small splitting of the rotational levels arising through 
lambda doubling and hyperfine splitting. This enables photons emitted in one infrared 
transition between sub levels of different rotational states in one OH molecule to to be 
absorbed by a transition between different sub-levels in the same pair of rotational rates. 
This tends to scramble the level populations within the upper or lower pairs of sub levels 
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Figure 2. Grey-scale shows the continuum emission of IC 443 at 1667 MHz (Jy per 76" x 51" 
beam). The beam size and orientation are indicated by the ellipse in the lower left corner. Con- 
tours show the velocity-integrated OH absorption from 0.2 to 1.6 Jy beam~^ kms~^ in steps 
of 0.2 Jy beam~^ kms~^. The three crosses mark positions of 1720 MHz OH maser emission 
(Hewitt 2009). The OH absorption clumps B, C, D, F, and G are labelled following the identifi- 
cation of shocked CO clumps by Huang, Dickman & Snell (1986). The clump of OH absorption 
labelled OHl A has not been previously identified. 



Table 2. SNR with iVoH > 10^'' cm"^ inferred from GET observations. 



SNR 


iVoH (dex) 


6049 MHz 


iVoH (dex) 




from GET 


searches^ 


from VLA 


W28 


17.14 ± 0.01 


PEA 




G16.7-f0.1 


16.9 ±0.6 


PEA 




3C391 


17.1 ± 0.4 


P E 


< 16.5 


W44 


17.14 ± 0.01 


P E 


< 17.1 


1C443 


17.11 ± 0.01 


P E 


~ 17.5 


CTB37A 


17.1 ± 0.4 


P A 




G349.7-H0.2 


17.04 ± 0.02 


P 




G357.7-0.1 


17.12 ±0.01 


P A 




G359.1-0.5 


17.13 ±0.01 


P A 





1 P=Parkes, E=Effelsberg, A=ATCA 



inside each rotational level and suppress population inversion in the 1720 and 6049 MHz 
transitions, quenching the masers. 

Line overlap is at present treated locally in the OH excitation code, i.e. the OH-bearing 
cloud is treated as a uniform slab with a velocity width that is everywhere the same. 
Its effects are illustrated in Fig. |4] The upper panel shows the maser optical depth as 
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RA (J2000) 

Figure 3. Grcy-scalo shows the continuum emission of the NE corner of 3C391 at 1667 MHz 
between -0.1 and 2.7 Jy per 96" x 56" beam. Absorption in the 1667MHz hue of OH is indicated 
by the contours, with levels at —1, —0.8, —0.6, —0.4, —0.2 Jy beam"'^ km s^^. This image is 
~7.44 arcmin on a side, approximately matching the GBT beam area of the 3C 391 spectrum 
from Hewitt et al. (2008). 



a function of OH column density in a uniform slab when overlap is neglected. In this 

case, line optical depths (ic. photon escape probabilities) scale as Nqk/^v so the curves 
simply shift to the right as Av is increased. The right hand column shows the effect of 
including local line overlap: once the line width exceeds 0.5 kms~^, line overlap starts to 
shuffle the populations within rotational levels and the inversions are suppressed. 

Observed 1720 MHz maser lines have Av ~ lkms~^, at which point local line overlap 
would suppress the masers, so this treatment overestimates the effect of line overlap. 
In reality the OH column docs not have a single, uniform, velocity profile throughout 
- instead, the column of OH has local line widths that are much smaller than the net 
velocity width created by the systematic velocity gradient along the line of sight. This 
situation requires a "non-local" treatment of line overlap effects, in which photons emit- 
ted by an infrared transition are absorbed in a neighbouring transition at a different 
layer in the slab (e.g. Docl, Gray & Field 1990). This form of overlap is less effective, 
potentially enabling the 1720MHz line to remain inverted for larger (~ lkms~^) line 
widths; the 6049 MHz line however, may then only be weakly inverted at column densi- 
ties ^ lO^^'cm"^. Resolution of this will have to await calculations including non-local 
overlap. 
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Figure 4. 1720 MHz (red) and 6049 MHz (blue) maser optical depth as a function of OH 
column density for a uniform slab with n(H2) — 10*cm~^ and T = 50K as the total (thermal 
+ micro-turbulent) FWHM line width is increased: O.lkms"^ (solid), 0.5kms~^ (long-dashed), 
0.6kms~^ (short- dashed), 0.7kms~^ (dotted). The upper and lower panels show the results 
obtained neglecting or including local line overlap, respectively (see text). 

4. Discussion 

The lack of detected OH(6049MHz) masers in the SNR searches conducted to date 
may be attributed to a combination if two factors. First, the OH column densities inferred 
from GET observations may be systematically overestimated because of the tendency for 
OH absorption to overlay regions where the background continuum is brightest. If the 
continuum has structure within the beam the optical depth (and so the OH column) will 
tend to be overestimated, as appears to be the case in 3C 391 and W44. However, IC 443 
clump G appears to have a column density 3 x 10^'' cm~^, more than twice that estimated 
from the single-dish observations. Therefore several additional remnants listed in Table 1 
may still have OH column densities exceeding the ~ 10^^ cm~^ threshold suggested for 
6049 MHz inversion. Second, non-local line overlap might suppress 6049 MHz inversion 
more effectively than the 1720 MHz transition, or simply push the threshold OH column 
density required for significant maser amplification to columns exceeding 3 x 10^^ cm~^. 

Acknowledgements 

We thank Jack Hewitt for providing his 1720 MHz VLA data. 



8 



Mark Wardle & Korinne McDonnell 



References 

Arikawa, Y., Tatematsu, K., Sekimoto, Y., & Takahashi, T. 1999, PASJ, 51, L7 

Brogan, C. L., Frail, D. A., Goss, W. M., & Troland, T. H. 2000, ApJ, 537, 875 

Claussen, M. J., Frail, D. A., Goss, W. M., & Gaumc, R. A. 1997, ApJ, 489, 143 

Claussen, M. J., Goss, W. M., Frail, D. A., & Desai, K. 1999, ApJ, 522, 349 

Cox, D. P., Shelton, R. L., Maciejewski, W., Smith, R. K., Plewa, T., Pawl, A., & Rozyczka, M. 

1999, ApJ, 524, 179 
DcNoycr, L. K. 1979a, ApJ (Letters), 228, L41 
DcNoycr, L. K. 1979b, ApJ (Letters), 232, L165 
Doel, R. C., Gray, M. D., & Field, D. 1990, MNRAS, 244, 504 
Drains, B. T. 1980, ApJ, 241, 102 

Draine, B. T., Roberge, W. G., & Dalgarno, A. 1983, ApJ, 264, 485-507 
Elitzur, M. 1976, ApJ, 203, 124 

Fish, V. L., Sjouwerman, L. O., & Pihlstrom, Y. M. 2007, ApJ (Letters), 670, L117 
Frail, D. A., Goss, W. M., & Slysh, V. I. 1994, ApJ, 424, Llll 

Frail, D. A., Goss, W. M., Reynoso, E. M., Giacani, E. B., Green, A. J., & Otrupcek, R. 1996, 
AJ, 111, 1651 

Goss, W. M. & Robinson, B. J. 1968, Astrophys. Lett., 2, 81 
Green, A. J., Frail, D. A., Goss, W. M., & Otrupcek, R. 1997,^J, 114, 2058 
Green, A. J., Wardle, M, & Lazcndic, J. S. 2000, lAU XXIV General Assembly, JDl, abstract. 
Hewitt, J. W., Yusef-Zadeh, F., Wardle, M., Roberts, D. A., & Kassim, N. E. 2006, ApJ, 652, 
1288 

Hewitt, J. W., Yusef-Zadch, F., & Wardle, M. 2008, ApJ, 683, 189 
Hewitt, J. W., Yusef-Zadch, F., & Wardle, M. 2009, ApJ (Letters), 706, L270 
Hewitt, J. W., 2009, PhD thesis. Northwestern University 
Huang, Y.-L., Dickman, R. L., & Snell, R. L. 1986, ApJ (Letters), 302, L63 
Hoffman, I. M., Goss, W. M., Brogan, C. L., Claussen, M. J., & Richards, A. M. S. 2003, ApJ, 
583, 272 

HoUonbach, D. J. & McKcc, C. F. 1989, ApJ, 342, 306 

Koralesky, B., Frail, D. A., Goss, W. M., Claussen, M. J., & Green, A. J. 1998, AJ, 116, 1323 
Kaufman, M. J. & Neufeld, D. A. 1996, ApJ, 456, 250 

Lazendic, J. S., Wardle, M., Whiteoak, J. B., Burton, M. G., & Green, A. J. 2010, MNRAS, 

409, 371 

Lockott, P., Gauthier, E., & Elitzur, M. 1999, ApJ, 511, 235 

Long, K. S., Blair, W. P., Matsui, Y., & White, R. L. 1991, ApJ, 373, 567 

McDonnell, K. E., Wardle, M., & Vaughan, A. E. 2008, MNRAS, 390, 49 

McDonnell, K. E., 2011, PhD thesis, Macquarie University 

Mullan, D. J. 1971, MNRAS, 153, 145 

Pavlakis, K. G. & Kylafis, N. D. 2000, ApJ, 534, 770 

Reach, W. T. & Rho, J. 1999, ApJ, 511, 836 

Rho, J. & Petre, R. 1998, ApJ (Letters), 503, L167 

Sjouwerman, L. O. & Pihlstrom, Y. M. 2008, ApJ, 681, 1287 

Wardle, M. 1999, ApJ (Letters), 525, LlOl 

Wardle, M. & Yusef-Zadeh, F. 2002, Science, 296, 2350 

White, R. L. & Long, K. S. 1991, ApJ, 373, 543 

Wootten, A. 1981, ApJ, 245, 105 

Yusef-Zadeh, F., Roberts, D. A., Goss, W. M., Frail, D. A., & Green, A. J. 1996, ApJ (Letters), 
466, L25 

Yusef-Zadeh, F., Roberts, D. A., Goss, W. M., Frail, D. A., & Green, A. J. 1999a, ApJ, 512, 230 
Yusef-Zadeh, F., Goss, W. M., Roberts, D. A., Robinson, B., & Frail, D. A. 1999b, ApJ, 527, 
172 

Yusef-Zadeh, F., Wardle, M., & Roberts, D. A. 2003, ApJ, 583, 267 
Yusef-Zadeh, F., Wardle, M., Rho, J., & Sakano, M. 2003, ApJ, 585, 319 



